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hydrogen  ion  due  to  lower-hybrid  heating  in  the  suprauroral  region 
based  on  a  mean  particle  calculation.  Initial  energy  of  the  ion  was 
taken  as  1  eV.  The  value  |Ey,,|2/G)LH  was  assumed  to  vary  linearly  from 
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and  field-aligned  directions.  The  frequencies  along  the  three  rays 
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8.  Typical  calculated  energetics.  The  lower  panel  gives  the  electrostatic  20 
energy  density  as  a  fimcdon  of  time,  while  die  top  panel  gives  the  ion 
energy  density.  The  solid  curves  are  calculated  based  on  die  theoretical 
model  of  RCJ.  The  dot-dashed  curves  give  the  corresponding  results 

when  mode-coupling  is  neglected.  The  dashed  curve  in  the  bottom 
panel  gives  the  result  from  an  equivalent  PIC  simulation. 

9.  Typical  calculated  ion-velocity  distributions.  The  solid  curve  is  a  21 
snapshot  of  the  velocity  distribution  at  a  time  when  mode-coupling 
process  has  taken  place.  The  da^ed  curve  is  obtained  from  a 
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based  on  the  RCJ  dieory.  Shading  indicates  the  strengdi  of  the  electric 
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iield  amplitude.  Notice  how  the  initial  nearly  uniform  wave  amplitude 
intensity  evolves  into  intense  localized  solitary  structures. 

12.  Typical  observed  spikelets  from  the  MARIE  rocket.  Top  panel  is  a  22 
frequency  spectrum  and  bottom  panel  shows  the  electric  field 
wavefonn.  (Courtesy  of  J.  LaBelle  and  P.M.  Kintner;  from  Ref  8.) 

13.  An  example  of  a  typical  spikelet  from  the  MARIE  rocket  at  an  altitude  2  3 
of  593  km.  The  upper  panel  shows  50  ms  of  electric  field  data  and 
bottom  panel  is  an  expanded  view  of  5  ms  of  data  centered  at  t  = 
249.689.  (Courtesy  of  J.  LaBelle  and  P.M.  Kintner;  from  Ref  8.) 

14.  Two  seconds  of  TOPAZ  3  data.  Upper  panels  give  the  count  rate  2  3 
summed  over  all  pitch  angles  and  energy  sweep  data.  Second  panel 
from  the  top  shows  the  ion  count  rate  as  a  fimction  of  pitch  angle. 

Third  panel  from  the  top  shows  the  electric  field  power  spectrum  as  a 

gray  scale.  Bottom  panel  shows  the  current  drawn  by  a  fixed  bias 
(+5V)  Langmuir  probe.  (Courtesy  of  P.M.  Kintner  and  J.  vago;  from 
Ref  8.) 

15.  Another  two-seconds  of  data  for  the  TOPAZ  3  rocket.  Same  format  as  24 
Figure  14.  (Courtesy  of  R.L.  Amoldy  and  K.  Lynch;  from  Ref  11.) 

16.  Phase  space  plots  for  the  electrons,  oxygen,  and  hydrogen  ions.  Top  24 
panels  are  for  a  time  when  transversely  accelerated  ions  were  not 
measured.  Bottom  panels  are  obtained  when  transversely  accelerated 

ions  were  detected.  The  contours  are  for  constant  log  f  and  the  heavy 
lines  space  by  one  unit  of  log  f  "r  has  units  of  sec^km"^.  (Courtesy  of 
R.L.  Amoldy  and  K.  Lynch;  from  Ref.  1 1.) 

17.  Velocity  phase  space  density  of  900-120°  transversely  accelerated  ions  25 
from  the  MARIE  rocket.  Dashed  curves  are  Maxwellian  fits. 
(Courtesy  of  A.W.  Yau  and  B.A.  Whalen;  from  Ref.  9.) 

18.  Monte  Carlo  simulation  of  an  ion  conic  for  parameters  consistent  with  25 
that  of  the  MARIE  flight.  The  normalization  of  phase  space  density  is 
arbitrary. 

19.  Counterstreaming  electrons  observed  by  the  DE2  satellite.  (Courtesy  of  26 
C.S.  Lin,  from  ref.  7.) 

20.  Electron  velocity  distribution  in  a  two-dimensional  simulation.  The  26 
dash-dot  curves  give  the  initial  parallel  velocity  distribution.  Dashed 

and  solid  curves  are  those  at  two  times  near  the  end  of  the  simulation 
run. 
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SYNOPSIS  OF  ACTIVITIES 


During  the  past  year,  we  have  made  substantive  progress  toward  the 
goals  set  forth  in  our  research  proposal:  "Mesoscale  Ionospheric 
Phenomena."  In  the  research  area,  we  have  succeeded  in  developing  the 
basic  understanding  of  the  mesoscale  processes  of  lower  hybrid  collapse, 
caviton  turbulence  and  charged  particle  energization  in  the  topside 
ionosphere  and  magnetosphere.  We  have  also  made  considerable  advance 
in  die  study  of  the  polar  wind  and  ion  cyclotron  resonance  heaing  of  central 
plasma  sheet  heavy  ions.  A  number  of  these  theoretical  findings  have 
recently  been  confirmed  by  experimental  data  collected  by  polar-orbiting 
satellites  as  well  as  high  altitude  sounding  rockets.  As  a  gauge  of 
recognition  of  the  group's  quality  of  research  by  its  peers,  our  group 
presented  a  total  of  16  invited  and  review  lectures  during  die  year  at  various 
national  and  international  conferences  and  research  institutions. 

In  January  of  1992,  we  organized  an  annual  Symposium  on  the 
"Physics  of  Space  Plasmas".  Professor  Eugene  Parker  of  die  University  of 
Chicago  was  the  Awardee  of  the  1992  Alfven  Lectureship  and  delivered  an 
Opening  Lecture  on  the  subject  matter  of:  "Spontaneous  Discontinuities  and 
Stellar  X-Rays."  We  also  organized  a  1992  Cambridge  Workshop  on 
Theoretical  Geoplasma  Research.  The  theme  of  the  workshop  was 
"Controversial  Issues  and  New  Frontier  Research  in  Geoplasmas."  This 
activity  was  participated  by  over  125  research  scientists  and  graduate 
students  fi'om  all  comers  of  the  world.  Ehiring  the  year,  we  also 
participated  in  the  scientific  sessions  at  the  Spring  and  Fall  AGU  Meetings, 
the  Plasma  Physics  Meeting  of  the  American  Physical  Society,  the  Third 
Huntsville  Workshop  on  Magnetosphere/Ionosphere  Plasma  Models,  the 
Western  Pacific  Geophysics  Meeting,  the  Nineteenth  IEEE  International 
Conference  on  Plasma  Science,  the  Chapman  Conference  on  Micro-Meso 
Scale  Phenomena  in  Space  Plasmas,  and  die  28th  Annual  Meeting  of  the 
Society  of  Engineering  Science. 

In  this  report,  we  provide  a  comprehensive  review  of  our  contribution 
toward  the  understanding  of  lower  hybrid  ion  heating  processes  in  Earth's 
ionosphere  and  magnetosphere.  This  review  is  prepared  by  Tom  Chang  for 
the  fordicoming  Special  Issue  of  Invited  Lectures  of  Physics  of  Fluids  B, 
Vol.  5. 
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1.  INTRODUCTION 


In  this  review,  we  consider  an  interesting  plasma  physics  phenomenon  that 
occurs  rather  frequently  within  the  Earth's  magnetosphere.  The  Earth's 
magiietosphere  is  a  cavity  carved  out  from  the  solar  wind  by  the  Earth’s  magnetic 
dipole  (Figure  1).  The  plasma  density  of  the  magnetosphere  is  quite  low  (of  the 
order  of  one  particle  per  cubic  centimeter)  and  most  of  the  particle  population  of 
the  magnetosphere  resi^’es  in  a  region  called  the  "Plasma  Sheet."  The  plasma 
sheet  is  quite  dynamic,  expanding  and  contracting  aperiodically  in  time. 
Frequently,  field-aligned  energetic  (keV)  electrons  have  been  detected  streaming 
toward  the  Earth's  ionosphere  near  the  boundary  layer  region  of  the  plasma  sheet 
(Figure  2).  In  addition  to  being  responsible  for  the  production  of  the  visible 
discrete  auroral  display  in  the  E-region  of  the  ionosphere,  these  energized 
electrons  are  also  the  culprit  for  a  suite  of  plasma  waves  that  are  observed  along 
the  auroral  field  lines.  In  particular,  intense  electrostatic  lower  hybrid  turbulence 
has  been  detected  ^  near  the  topside  ionosphere  and  the  suprauroral  region  where 
the  plasma  beta  is  very  low  (Figure  3). 

The  particle  populations  in  the  suprauroral  region  are  generally  quite 
anisotropic.  For  example,  ion  distributions  strongly  peaked  in  pitch  angle  are 
routinely  observed^  (Figure  4).  The  conic  shape  of  these  distributions  indicates 
some  form  of  heatmg  transverse  to  the  ambient  magnetic  field.  The  generally 
accepted  scenario  for  such  transverse  acceleration  is  some  sort  of  wave-particle 
interaction.  In  this  picture  (Figure  5),  the  ions  are  energized  perpendicularly  to  the 
geomagnetic  field  lines  by  energy-carrying  plasma  waves.  One  can  then  account 
for  the  conic  form  of  the  distribution  by  realizing  that  the  magnetic  field  strength 
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decreases  with  altitude.  Thus,  the  adiabatic  motion  of  the  ions  drifting  to  higher 
altitudes  transforms  the  heated  distribution  into  one  that  is  more  fteld-aligned,  i.e., 
a  conic^. 

It  was  suggested  by  Chang  and  Coppi^  in  1981  that  the  intense  lower 
hybrid  waves  that  are  detected  along  the  discrete  auroral  field  lines  could  be  the 
prime  cause  for  the  ion  acceleration  process.  It  is  known  that  electron  beams  such 
as  the  precipitating  energetic  auroral  electrons  can  provide  the  free  energy  to 
generate  waves  in  the  whistler  range  of  fi’equencies  (VLF  waves  for 
suprauroral/auroral  parameters)  including  ftie  lower  hybrid  waves^.  Since  the 
observed  lower  hybrid  waves  were  generally  broad  band,  a  quasilinear  diffusion 
operator  was  used  to  estimate  the  average  energy  transfer  from  a  steady  state  of 
waves  populating  a  portion  of  the  field  line.  It  was  found  that  typically  1-eV  ions 
could  be  raised  in  energy  to  tens  or  hundreds  of  eV  and  beyond  by  lower  hybrid 
waves  of  moderate  intensity  along  the  discrete  auroral  field  lines  where  electron 
beams  were  detected  (Figure  6).  As  the  heated  ions  drift  upward  along  the  field 
lines  due  to  the  "mirror"  geometry  of  the  Earth's  magnetic  field,  they  eventually 
leave  the  primary  heating  region.  Thus,  the  heating  process  was  found  to  be  self- 
limiting. 

However,  we  are  confronted  with  the  dilemma  that  the  lower  hybrid  waves 
initially  excited  by  the  linear  instability  of  the  energetic  auroral  electrons  are  of 
high  phase  velocity,  they  would  resonantly  interact  only  with  energetic  ions, 
which  are  few  in  number  in  the  topside  ionosphere.  In  this  review,  we  shall 
describe  a  scenario  (first  introduced  by  Retterer,  Chang,  and  Jasperse  (RCJ)^  in 
1986)  utilizing  the  convective  nature  of  the  waves  and  ensuing  nonlinear  mode- 
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coupling  processes  that  can  result  in  the  collapse  of  the  lower  hybrid  modes  to 
shorter  wavelengths  to  achieve  the  resonance  matching  with  the  tail  region  of  the 
cold  ionospheric  ion  distribution.  The  theory  also  predicts  the  simultaneous 
occurrence  of  counterstreaming  electrons  that  are  commonly  observed^  in 
conjunction  with  the  ion  conics  in  the  suprauroral  region.  These  theoretical 
predictions,  including  the  eventual  generation  of  lower  hybrid  caviton  turbulence 
due  to  the  phenomena  of  collapse  throi^  modulational  instability  have  been 
confirmed  by  recent  high-altitude  sounding  rocket  experiments*"!^,  utilizing 
iimovative  high-time  resolution  wave  and  particle  detectors. 

2.  THE  THEORY  OF  RETTERER,  CHANG,  AND  JASPERSE  (RCJ) 

The  commonly  observed  geometry  of  discrete  aurorae  in  the  high  latitude 
ionosphere  is  that  the  auroral  arcs  are  sufficiently  extended  in  the  longitudinal 
direction  but  narrowly  confined  in  the  latitudinal  direction.  This  indicates  that  the 
regime  containing  the  energetic  precipitating  electrons  above  a  discrete  auroral  arc 
(the  suprauroral/auroral  region)  must  also  be  similarly  confined  in  the  latitudinal 
direction.  The  density  within  the  suprauroral  region  is  generally  less  than  that  of 
its  immediate  neighborhoods  with  appreciable  perpendicular  (outward)  and 
parallel  (downward)  density  gradients.  The  propagation  characteristics  of  the 
relevant  plasma  waves  that  are  excited  by  the  energetic  electron  beam  at  high 
altitudes  depend  crucially  on  the  inhomogeneous  density  structure  of  the  depleted 
suprauroral  region.  It  has  been  shown  by  Maggs^"!,  Maggs  and  Lotko^^,  and 
Retterer  et  al.!^,  using  ray-tracing  calculations,  that  for  typical  plasma  parameters 
in  the  suprauroral  region,  convective  modes  at  frequencies  much  higher  than  the 
lower  hybrid  frequency  will  generally  prope^ate  out  of  the  region  in  the  latitudinal 
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direction  shortly  after  thev'  are  generated,  leaving  only  those  modes  close  to  the 
lower  hybrid  frequency,  which  have  group  velocities  nearly  parallel  to  the 
magnetic  field,  to  continue  to  propagate  and  grow  along  die  field-aligned 
direction.  These  lower  hybrid  modes  can  stay  withm  the  suprauroral  region  and 
grow  to  very  large  amplitudes  leading  to  die  nonlinear  coupling  of  modes  and 
collapse  of  waves  to  shorter  wavelengdis  near  the  topside  ionosphere  (Figure  7) 

It  has  been  demonstrated  by  Retterer.  Chang  and  Jasperse^  (RCJ)  that  such 
type  of  nonlinear  coupling  process  may  be  understood  in  terms  of  the  phenomenon 
of  modulational  instability  based  on  the  warm  fluid  approximation  for  the  ions  and 
electrons.  (Similar  ideas  have  been  considered  by  odier  authors  in  other  plasma 
applications  and  in  the  laboratory  Assuming  that  the  high  frequency 
(lower  hybrid)  fluctuations  arc  modulated  by  the  low  frequency  (acoustic) 
components,  the  resulting  equation  governing  the  amplitude  of  die  lower  hybrid 
electric  field  for  propagations  restricted  at  an  angle  characterized  by  die  square 
root  of  the  mass  ratio  (k*B  » (m^mj)^^)  is  a  cubic  nonlinear  Schrbdinger  equation 
(NLSE).  (The  assumption  of  one-dimensional  propagations  is  a  physically 
realistic  approximation.  Due  to  the  narrowness  of  the  suprauroral  region  in  the 
latitudinal  direction,  other  nonlinear  mode-coupling  processes  diat  are  neglected  in 
this  approximation  entail  die  ExB  drift  of  the  electrons  and  dierefore  can  not 
support  the  continued  growth  of  the  ensuing  instabilities.  The  choice  of  the 
particular  angle  given  above  facilitates  the  growth  of  waves  near  the  lower  hybrid 
frequency.) 

In  addition  to  the  nonlinear  coupling  process  of  wave-wave  interactions, 
concomitant  wave-particle  interactions  control  the  evolution  of  the  lower  hybrid 
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turbulence  in  oSe  tc^ide  ionosphere  RCJ*  suggested  that  the  dominant  wave- 
particle  plasma  interaction  process  could  be  modeled  bv  the  consideration  of 
instantaneous  Landau  damping  of  the  waves  by  the  various  particle  species  (i  e  . 
the  beam  electrons,  the  background  ambient  ions  and  electrons)  which  at  the  same 
time  are  scattered  by  the  lower  hybrid  waves  through  resonant  interactions  The 
resulting  coupled  equations  characterizing  both  the  nonlinear  wave- wave  and 
wave-particle  interacticMis  given  below  have  been  solved  numerically*  for  typical 
suprauroral  conditions  near  the  topside  ionosphere 

idEJdt  =  Dk}E, -C'^E,.e:.E,,,..,  +/>,(*)£»  (1) 

af,(y)idt  =  diav[D,iv)^jav]  (2) 

where  £*  is  die  Fourier  component  of  the  lower  hybrid  electric  field  fluctuations, 
C  and  D  are  die  coupling  and  diffusion  coefficients  expressible  in  terms  of  the 
plasma  parameters,  is  the  instantaneous  Laudau  damping  rate,  and  are  the 
distribution  functions.  Because  the  average  intensities  of  the  waves  are 
sufficiendy  low  and  die  spectrum  is  broad-band,  RCJ  used  die  quasilinear 
diffiision  operator  for  die  scattering  of  particles  by  die  waves.  This  approximation 
is  quite  adequate  except  during  die  initial  transitional  regime  firom  linear  to  die 
saturation  state.  It  is  to  be  noted  that  using  the  quasilinear  diffusion  operator  is  not 
the  same  as  assuming  quasilinear  plateauing.  The  turbulence  is  strong  and  is 
governed  by  Eq.(l).  Quasilinear  theory  would  have  entailed  a  linear  growth 
equation  for  E]. 
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The  bottom  panel  of  (Figure  8)  gives  the  calculated  evolution  of  the  lower 
hybrid  wave  energy  density  as  a  function  of  elapsed  time.  It  is  seen  diat  the  wave 
energy  grows  linearly  during  the  initial  stage  of  excitation  and  dien  saturates  due 
to  mmlinear  mode-coupling  processes.  The  result  compares  favorably  widi  that 
obtained  from  an  equivalent  PIC  simulation.  The  top  panel  of  (Figure  8)  gives  the 
time  evolution  of  the  ion  energy  density.  It  is  noted  that  the  ions  are  heated  only 
after  nonlinear  mode  coupling  has  set  in.  At  later  stages  of  die  evolution  when  the 
waves  are  saturated,  die  ions  are  seen  to  continue  to  gain  dieir  energy.  The  lower 
hybrid  waves  now  simply  play  die  role  of  a  conduit,  transferring  energy  from  the 
electron  beam  to  the  ions.  If  mode-coupling  processes  are  not  allowed  to  take 
place  (quasilinear  dieory),  dien  die  waves  will  continue  to  grow  until  quasilinear 
plateauing  of  die  beam  electrons  sets  in  and  the  ions  cannot  be  heated.  Figure  9 
compares  the  ion  velocity  distribution  at  die  end  of  a  PIC  simulation  with  the 
dieoredcal  prediction  based  on  the  hybrid  kinetic/NLSE  model  of  RCJ.  Note  that 
the  tail  regions  of  an  initial  Maxwellian  ion  distribution  have  been  heated  by  die 
nonlinear  coupling  process. 

Figure  10  presents  the  evolution  of  the  wave  spectrum  in  a  one-dimensional 
PIC  simulation.  It  is  seen  that  the  wave  energy  cascades  continually  from  the  k- 
values  of  the  linearly  unstable  modes  toward  larger  k-values  indicating  the 
collapsing  process  toward  shorter  wavelength  lower  hybrid  waves.  Figure  11 
depicts  the  evolution  of  the  amplitude  of  the  electric  field  in  real  space  and  time 
based  on  die  RCJ  model^O.  At  the  beginning,  the  amplitude  of  die  electric  field 
grows  uniformly.  Eventually,  when  wave  collapse  sets  in,  the  electric  field 
intensity  begins  to  concentrate  aperiodically  and  in  narrow  spatial  regions.  These 
localized  structures  (cavitons,  because  they  correspond  to  sharp  localized  density 
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depletions)  are  seen  to  propagate  at  speeds  much  less  than  the  ambient  ion  thermal 
velocity  and  should  appear  as  stationary  structures  if  detected  by  in  situ  electric 
field  measurements.  The  onset  of  caviton  turbulence  b  expected  from  the 
model  equations  since  it  b  known  that  the  cubic  NLSE  can  support  such  type 
of  solitary  structures'^ 

3.  CONFIRMATION  OF  THE  RCJ  THEORY  BY  ROCKET 
OBSERVATIONS 

Recent  rocket  experiments*"*  3  (MARIE  and  TOPAZ  3)  carried  out  by 
Kintner,  Yau,  Labelle,  Amoldy,  and  their  co-workers  have  obtained  wave-particle 
data  in  die  topside  ionosphere,  that  are  in  remarkable  agreement  with  the 
predictions  of  die  RCJ  theory. 

For  example,  electric  field  data  from  both  the  MARIE  and  TOPAZ  3 
satellites  have  obtained  irrefutable  evidence  of  caviton  turbulence  in  the  topside 
ionosphere  widi  each  individual  caviton  exhibiting  the  characteristic  envelope 
solitary  structure  as  predicted  by  the  cubic  NLSE  and  die  RCJ  theory,  (Figures  12- 
13).  12-13.  The  intense  cavitons  are  clearly  seen  in  Figure  12  as  spiky  structures* 
with  high  intense  electric  field  intensity  (50-100  mV/m  for  the  MARIE  satellite). 
These  structures  have  been  christened  "Lower  Hybrid  Spikelets"  by  LaBelle  et  al.* 
A  typical  example  of  one  of  these  spikelets  is  displayed  in  expanded  time  scale  in 
Figure  13;  the  envelope  structure  with  characteristic  lower  hybrid  oscillations  is 
clearly  displayed*.  The  frequency  of  the  peak  of  the  spectral  density  of  one  of  the 
typical  spikelets  was  around  10  kHz  and  a  wavelength  of  approximately  20  m, 
giving  a  characteristic  phase  velocity  of  waves  at  200  km/s^^.  Similar  events  are 
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routinely  seen  by  the  TOPAZ  3  satellite*®’ *2  ^ith  peak  intensities  at  300-400 
mV/m  or  beyond  when  die  instrument  is  saturated. 

That  these  structures  are  cavitons  are  clearly  seen  in  Figures  14-15*®"** 
where  the  lower  panels  give  Langmuir  probe  measurements  indicating  highly 
rarefied  regions  corresponding  to  the  locations  of  die  electric  field  spikelets  (third 
panels  fix>m  the  top).  Simultaneous  with  die  observation  of  cavitons,  the  ions  are 
seen  accelerated  in  die  90  degree  pitch  angle  direction.  In  Figure  15,  one  of  die 
ion  events  (at  577. 1  seconds)  is  not  accompanied  by  the  simultaneous  observation 
of  a  caviton  structure.  The  interpretation  by  Amoldy  et  al.  *  *>*^  is  diat  the  ions  are 
not  exactly  at  90  degrees  and  are,  therefore,  probably  coming  from  a  region  below 
the  ^cecraft.  This  indicates  that  it  is  the  cavitons  that  accelerate  the  ions  radier 
than  die  ions  producing  the  waves;  for  odierwise,  cavitons  should  always  be 
present  when  transversely  accelerated  ions  are  seen. 

Figure  16  gives  typical  particle  distributions  obtained  by  the  TOPAZ  3 
satellite**  The  top  diree  panels  are  phase  space  plots  when  transverse  heating  is 
not  observed.  The  lower  panels  are  phase  plots  vdien  transverse  heating  is 
observed.  At  t  =  577.5  sec.,  the  oxygen  population  is  tailed  heated  widi  the  tail 
emanating  from  vj^  approximately  equal  to  10  km/s.  Simultaneously,  the  hydrogen 
population  is  also  heated.  For  lower  hybrid  heating,  the  hydrogen  population 
should  again  be  heated  from  vj^  =  10  km/s.  This  is  consistent  from  the  near  bulk 
heating  signature  of  the  phase  space  plot.  We  defer  the  discussion  of  electron 
signatures  in  a  later  section. 

Generally,  the  ions  will  receive  many  kicks  by  the  cavitons.  The  integral 
effect  on  the  ion  distribution  by  the  ensemble  of  cavitons  is  the  diffusion  of  ions  in 


8 


velocity  space,  which  can  be  estimated  from  the  ensemble  statistics  of  the 
numerical  solution  of  Eqs.  (1)  and  (2).  Alternatively,  the  perpendicular  velocity 
diffusion  coefficient  may  be  estimated  from  die  saturated  k-spectrum  such  as  those 
given  in  Figure  10.  Based  on  such  spectrum  forms.  Crew  and  Chang^^  suggested 
the  following  form  for  the  velocity  dependence  of  the  diffusion  coefficient; 

where  is  a  reference  value  of  die  electric  field  spectral  intensity,  v,,  is  the 
reference  phase  velocity  and  ^  is  a  parameter  which  depends  on  the  shape  of  the 
spectrum  curve. 

Recendy,  Retterer,  Chang  and  Jasperse^O  used  the  above  form  of  the  difiiision 
coefficient  to  estimate  the  energization  of  ions  by  the  caviton  turbulence  observed 
by  the  MARIE  satellite*'^.  The  values  of  and  are  chosen  as  (7  mV/m)2/kHz 
and  200  km/s,  respectively  based  on  estimates  from  typical  wave  measurements  on 
MARIE.  Using  Monte  Carlo  techniques  originally  developed  by  Retterer,  Chang 
and  Jasperse^^  in  1983  and  with  a  judicious  choice  of  initial  conditions,  it  was 
shown  that  the  simultaneously  observed  hydrogen  data^«20  may  5^  effectively 
reproduced,  provided  the  shape  parameter  p  is  chosen  to  be  unity  (Figures  17-18). 

4.  ELECTRON  SIGNATURES 

Turning  our  attention  to  the  electrons,  we  first  note  that  the  collapse  process 
characterized  by  Eqs.  (1)  will  generate  short  wavelength  lower  hybrid  waves  with 
phase  velocities  symmetrically  in  both  positive  and  negative  directions  along  the 
angle  characterized  by  k«B  =  Thus,  we  expect  the  ambient  electrons 
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to  be  energized  bod)  upward  and  downward  along  the  field  lines.  This  type  of 
counterstreaming  electrons  have  been  observed  in  the  magnetosphere  in 
conjunction  with  ion  conics^  (Figure  19).  In  addition,  of  course,  we  expect  to  see 
a  flattened  (N.B.,  not  quasilinear  plateaued)  component  of  the  beam  population. 
However,  unlike  die  ions,  the  electron  population  is  also  affected  by  VLF  or  lower 
hybrid  waves  at  angles  other  dian  diat  characterized  by  die  square  root  of  the  mass 
ratio.  Thus,  a  two-dimensional  simulation  or  analysis  is  needed  to  study  die  full 
effect  of  the  VLF  waves  on  the  electron  population.  Such  types  of  simulations 
have  been  performed  by  Retterer  et  al.  1^,25  (figure  20). 

The  calculated  electron  population  is  very  similar  to  that  is  displayed  in  the 
left  panels  of  Figure  15  of  the  TOPAZ  data..  We  expect  to  see  die  flattened 
electron  beam  population  nearly  everywhere  within  die  discrete  auroral  region  and 
not  solely  in  conjunction  with  the  observation  of  cavitons.  Same  holds  true  for  the 
bi-directional  energization  of  the  ambient  population.  Measurements  have  also 
detected  short-lived  (the  order  of  seconds)  dispersive  bursts  of  cold  electrons  anti¬ 
correlated  with  the  transversely  accelerated  ions.  Their  role  in  die  totality  of  the 
wave-pardcle  interaction  process  within  a  discrete  auroral  arc  is  not  clear.  The 
flattened  electron  population  in  Figure  15  also  has  a  pitch-angle  scattered 
component.  This  could  be  accomplished  by  anomalous  Doppler  resonance  as 
suggested  by  Kadomtsev  and  Pogutse^^  and  Parail  and  Pogutse27.  n  has  been 
si^ested  by  Omelchenko  et  al.28  that  such  type  of  fan-like  electron  distribution 
could  enhance  the  local  production  of  lower  hybrid  waves. 

5.  SUMMARY  AND  RECENT  ACTIVITIES 
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In  this  review,  we  have  described  an  extremely  interesting  plasma  physics 
phenomenon  in  the  magnetosphere/ionosphere;  namely,  transverse  ion  acceleration 
by  localized  intense  lower  hybrid  cavitons.  The  theoretical  research  was  initiated 
in  198  H  in  response  to  the  curious  observations  of  ion  conics  in  the  suprauroral 
region  of  discrete  arcs.  This  was  followed  by  a  daring  theoretical  prediction 
(RCJ)^  which,  by  noting  that  the  relevant  plasma  process  is  micro/meso  scale  in 
nature,  suggested  that  all  the  important  actions  must  initiate  from  the  topside 
ionosphere.  The  theoretical  predictions  were  verified  by  a  series  of  recent 
innovative  rocket  experiments*"^^  carried  out  in  die  source  region.  Definitive 
evidence  of  die  existence  of  intense  lower  hybrid  spikelets  (or  cavitons)  were 
detected.  Concrete  evidence  was  collected  indicating  that  it  is  these  cavitons  that 
initiate  die  transverse  energization  of  ionospheric  ions  in  a  discrete  aurora. 
Ensemble-averaged  diffusion  calculations*^’20  of  die  ions  in  this  region  gave 
further  credibility  of  die  feasibility  of  energizing  ionospheric  ions  to 
magnetospheric  energies  by  broad  band  lower  hybrid  waves. 

Since  the  reporting  of  the  in-situ  measurements  of  cavitons  and  transverse 
ion  acceleration,  a  number  of  new  developments  have  emerged.  For  example, 
higher  dimensional  dieories  and  simulations  have  been  performed^  ^’25 ,29-31  n  jg 
now  understood  that  when  the  extra  dimensions  are  taken  into  consideration,  the 
cavitons  are  actually  cigar-shaped  structures  of  finite  extension.  The  major  axis  of 
each  individual  caviton  is  aligned  with  die  magnetic  field  and  the  aspect  ratio  is 
related  to  the  square  root  of  the  mass  ratio.  The  absolute  dimensions  of  these 
creatures  depend  on  the  details  of  the  wave-particle  interactions  and  the  ambient 
parameters.  Electron  dynamics  are  understood  more  fully  in  higher 
dimensions*^*20, 29-31  Limited  study  has  also  been  made  for  multi-ion 


11 


plasmas  Some  of  these  recent  calculations  are  performed  with  free  energies 
odier  than  that  of  the  electron  beam,  for  example  the  fan  instability  or  ion  ring 
distributions28.29,31,32  Applications  have  been  made  to  die  particle  energization 
processes  in  other  phenomena  in  space  plasmas;  e.g.,  cusp  ions^^  and  solar  flares 
where  a  relativistic  description  is  necessary^^.  It  has  also  been  suggested  that 
conversion  to  short  wavelength  lower  hybrid  waves  could  be  accomplished  by 
linear  conversions^  and/or  scattering  of  waves  over  prescribed  plasma 
irregularitiesS4. 

We  end  this  review  by  providing  three  additional  comments.  The  first 
comment  is  that  it  is  not  always  necessary  to  have  lower  hybrid  collapse  before 
transverse  ion  energization  can  take  place.  The  important  thing  is  that  there  is  a 
resonance  matching  between  die  phase  velocity  of  die  waves  and  the  particle 
velocity  of  the  ions.  This  point  is  demonstrated  rather  vividly  by  the  two  recent 
studies  of  Ergun  et  alS5  and  Ganguli  et  al.S6. 

The  second  comment  is  that  transverse  ion  energization  can  be 
accomplished  by  other  plasma  processes,  such  as  electromagnetic37-42  anj 
electrostatic^^  ion  cyclotron  resonance  heating,  multiple-cyclotron  resonance 
heating^,  oblique  double  layers^^^  stochastic  heating  for  coherent  waves^^,  and 
nonresonant  heating^^.  In  space  applications,  these  alternative  ideas  have  been 
reviewed  recendy  by  Lysak^^,  Chang  et  al.^^,  and  Chang  and  Andre^®,  among 
others. 


Finally,  most  of  the  particle  energization  processes  and  wave-particle 
interaction  phenomena  are  mesoscale  in  nature.  Self-consistent  calculations  diat 
address  both  the  background  inhomogeneous  plasma  and  the  microscopic 
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interactions  are  ultimately  needed.  Some  calculations  aiming  in  that  direction 
have  been  made  recently36,5l. 
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